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Abstract 

This paper discusses a few unique effects of ultra-thin- 
body double-gate NMOSFET that are arising from the 
bandstructure of the thin film Si channel. The bandstruc- 
ture has been calculated using 10-orbital sp 3 d 5 s* tight- 
binding method. A number of intrinsic properties includ- 
ing band gap, density of states, intrinsic carrier concentra- 
tion and parabolic effective mass have been derived from 
the calculated bandstructure. The spatial distributions 
of intrinsic carrier concentration and < 100 > effective 
mass, arising from the wavefunction of different contribut- 
ing subbands are analyzed. A self-consistent solution of 
Poisson-Schrodinger coupled equation is obtained taking 
the full bandstructure into account, which is then applied 
to an insightful analysis of volume inversion. The spatial 
distribution of carriers over the channel of a DGFET has 
been calculated and its effects on effective mass and chan- 
nel capacitance are discussed. 

1 Introduction 

The interest in Ultra-Thin-Body (UTB) Double-Gate 
FET (DGFET) has grown in the recent past because of 
its superior properties compared to bulk MOSFET and 
is being considered as one of the future alternatives of 
present day bulk devices [TJ [2] . Apart from superior gate 
control from both top and bottom, the intrinsic quantum 
confinement provided by its unique geometric structure 
affects the characteristics of UTB DGFET [3]. The dif- 
ferent aspects of classical modeling of DGFETs have been 
discussed in [H [5] [5] . There are a number of reports on 
the quantum mechanical effects on DGFET, both ana- 
lytical as well as numerical [5J [71 [HI E]. To solve nu- 
merically, one has to solve coupled Poisson-Schrodinger 
equations self-consistently [10] . This is done either by as- 
suming some analytical E — k relationship, or by taking 
full bandstructure into account. There has been consid- 
erable amount of work on calculation of bandstructure of 



materials [HI QJ H3 El H3 \M HZ] which can be plugged 
into the self-consistent Poisson-Schrodinger equation nu- 
merically [15] . 

DGFET has a unique property of volume inversion 
which improves the transport characteristics enormously 
PQ , [3] • This can be explained with the help of quantum 
effects. Another important aspect is a substantial change 
in transport properties depending on the crystallographic 
orientation [1^] [201 HI] which again can be analyzed from 
detailed bandstructure calculation. Also, the total num- 
ber of intrinsic carriers reduces with the thinning of the 
channel material. This has an effect on the total gate 
capacitance of DGFET [H [H] . 

The aim of this paper is to focus on detailed analysis 
of some effects in UTB-DGFET which arise entirely be- 
cause of bandstructure of the channel material and are 
not very apparent. Only Silicon has been considered as 
the thin channel material in this work, but this can be 
easily extended to other channel materials as well. The 
full-band structure calculation that has been used here is 
based on sp 3 d 5 s* tight-binding method [HI [03 [TB] . The 
calculated bandstructure has then been used to predict 
some intrinsic transport properties of thin film Si includ- 
ing band gap, density of states, intrinsic carrier concen- 
tration and effective mass. From this analysis, a number 
of features are explained which are unique to ultra-thin 
film semiconductors. Following this, Poisson-Schrodinger 
coupled equation is solved self-consistently taking care of 
the full bandstructure. With the help of this, volume 
inversion phenomenon has been critically analyzed with 
physical insights. This in turn throws some light on the 
spatial distribution of carriers inside the DGFET chan- 
nel. Taking this into account, total channel capacitance 
and evolution of effective mass from source end to drain 
end along the channel have been analyzed. 

The rest of the paper is organized as follows: Sec. [^dis- 
cusses on the details of sp 3 d 5 s* tight-binding method of 
bandstructure calculation. The different intrinsic trans- 
port properties of ultra-thin film Silicon have been dis- 
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cussed in sec. [3] Poisson-Schrodinger coupled equa- 
tion has been solved self-consistently and different related 
analyses have been performed in sec. [4j Finally the paper 
is concluded in sec. |H 



2 Bandstructure Calculation 

Tight-binding method of bandstructure calculation has 
been extensively studied by many researchers (TTJ [T2l H31 
HH HS1 HS1 H7]. In this work, a 10-orbital sp 3 d 5 s* tight- 
binding method O H3 [TB] has been used to find band- 
structure of the thin film Si, being used as channel ma- 
terial. Only the onsite energies and two-center overlap 
integrals of nearest neighbors have been taken into ac- 
count. Spin orbit interaction has been neglected, and 
thus each k point in the Brillouin zone is assumed to be 
degenerate with two spin states. Infinite crystal period- 
icity has been assumed along channel length and width 
directions and thus Bloch's theorem is assumed to hold 
good in those directions. However, along the thickness of 
the channel, the crystal is truncated to a few monolayers, 
thus crystal periodicity can not be assumed in this direc- 
tion. Suppose, the thickness contains N atomic mono- 
layers. Then, the truncated crystal can be formed by 
taking a basis of N atoms along the thickness direction 
and spanning them over the whole 2-D space. Fig. [I] 
shows a 7 monolayer thick channel with the basis atoms 
shown as black dots. The channel region can be formed 
by spanning the basis atoms along x and y. The tight- 
binding fitting parameters for Si, used in this work, have 
been taken from [15l [16] ■ An N monolayers thick film will 
produce a 107V x 107V tight-binding Hamiltonian [Tg]. To 
get rid of the huge number of surface states (whose energy 
eigen values often fall inside the semiconductor band gap) 
caused from dangling bonds, it has been assumed that the 
surfaces are completely passivated by Hydrogen. This has 
been achieved by artificially increasing the onsite energies 
of s and p orbitals of the surface atoms, as described in 

The assumption of this method is that the electronic 
wave function is strictly guided in the x — y plane. Thus, 
the Brillouin zone will comprise of a 2-D k space, as op- 
posed to a 3-D one in bulk case. k z has been assumed to 
be zero throughout this paper. The whole 2-D Brillouin 
zone has been discretized using step size of 0.05 x ^ for 
both k x and k y where a is lattice constant (=5.43A for 
Si). In this paper, the film thickness has been referenced 
to the number of monolayers (AL) in the film. An N 
AL thick Si film translates to a thickness of a(N — l)/4. 
Fig. [2] shows energy dispersion plot of a 17-monolayer 
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Figure 1: A 7-monolayer thick film with basis atoms 
(black dots). The basis atoms can be spanned in whole 
2 — D along x and y to construct the thin film. 



thick (~ 2.17nm) Si film over the whole 2 — D Brillouin 
zone. Only the top most valence subband and bottom 
most conduction subband have been included for clarity. 
Throughout this paper, the valleys occurring at T point 
and at ~ 0.8^ along X direction are termed as T valley 
and X valley respectively. 




Figure 2: E — k relationship of top most valence subband 
and bottom most conduction subband over the whole 2 — 
D Brillouin zone of a 17 monolayer thick Si film. The 
conduction band minimum occurs at V point. The X 
valley is 4 fold degenerate. 



2 



Intrinsic Properties 
Film Si 



of Thin 



In this section, the variations of different intrinsic electri- 
cal properties of thin film as a function of film thickness 
have been derived from the bandstructure calculation, de- 
scribed in the previous section. 

3.1 Bandgap and Density of States 

It has been well established in literature, by both theory 
and experiments, that in the nano-scale, bandgap of semi- 
conductors is a function of size of the material. As size 
reduces, the bandgap of the material increases. Fig. [3] 
shows how the T gap and X gap of a Si film are changing 
as a function of the film thickness. One should note that, 
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Figure 4: 2-D DOS as a function of electronic energy for 
different thickness values of the Si film. The energy space 
is discretized by steps of 0.2 eV. 



ent thickness values. In the calculation, the energy has 
been discretized in steps of 0.2eV. Just above the cut-off 
energy (conduction band minimum), only T valley con- 
tributes. However, as energy increases, other regions of 
Brillouin zone also start contributing. As expected, den- 
sity of states for thicker film is larger. 



Figure 3: Calculated T and X gap of Silicon thin film as a 
function of film thickness. For sufficiently small thickness, 
direct T gap is quite larger compared to the next gap 
occurring at X valley. 

for a sufficiently thin film, as opposed to the bulk case, the 
conduction band minimum occurs at direct T point, and 
not in the X direction. Thus the electrons will first popu- 
late the T valley and hence, one can expect to see drastic 
change in transport properties for a thin film Si channel 
compared to bulk. As the film thickness increases, the en- 
ergy difference AErx between T and X valleys decreases, 
and electron will start populating both the valleys. Fi- 
nally, at sufficiently large film thickness, at the bulk limit, 
X valley is of less energy compared to T valley, and thus, 
electrons will start populating only X valley. In Fig. |4j 
the 2-D density of states has been plotted as a function 
of electron energy in the conduction band, for 4 diffcr- 



3.2 Intrinsic Carrier Concentration 

As it has been discussed already that bandgap increases 
as the size goes from bulk to nano-scale, one expects lower 
intrinsic carrier concentration as the film thickness re- 
duces. The per unit area intrinsic electron concentration 
at temperature T is given by 



n A 



EE 2 M 



(i) 



where the first sum is over different subband indices j of 
conduction band and the second sum is over all k points 
in the first Brillouin zone. Ej represents the energy eigen- 
value at k th point of j th subband index. The Fermi-Dirac 
probability j{E^) is given by 



= 



1 
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ks is the Boltzmann constant and /i is the chemical po- 
tential. Fig. [5] shows that the intrinsic carrier concen- 
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Figure 5: Intrinsic carrier concentration per unit area (i.e. 
total number of carriers contained in the film with unity 
area) as a function of Si film thickness. 
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Figure 6: Percentage contribution to per unit area intrin- 
sic carrier concentration from different Si subbands lying 
in r and X valleys. For very small thickness only Ti 
and T 2 contribute, but at larger thickness, electrons start 
populating other valleys as well. 



tration per unit area increases with film thickness. How- 
ever, apart from the reduction in carrier concentration, 
another important observation is that the carriers have a 
distribution along the film thickness, which peaks at the 
center of the film. This is due to the spatial distribution 
of the wave function of the electronic states contributing 
to the carrier concentration. If the film of thickness t has 
N monolayers, then the film can be assumed to be dis- 
cretized by N points, at each of which, per unit volume 
carrier density is given by 

A*) = J ££2/(4)|^(z)| 2 (3) 

j k 

where ipj(z) is the wave function of the electronic state 
(j, k) at z. Fig. [6] plots the fractional contribution of dif- 
ferent subbands to total electron concentration for a thin 
film of Si. r, and Xj represent the i th subband of V valley 
and j th subband of X valley, respectively. It is clear that, 
for very small thickness, only T\ and T2 subbands con- 
tribute, but as thickness increases, other subbands also 
start contributing. In Fig. [7]and[8j the spatial distribu- 
tion of intrinsic carrier concentration, coming from differ- 
ent subbands, have been shown, for 9 and 33 monolayer 
thick Si film, respectively. Since for a 9 monolayer thick 
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Figure 7: Carrier distribution along channel thickness of 
different subbands lying in T and X valleys for a 9 atomic 
layer thick Si film. 
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Figure 8: Carrier distribution along channel thickness of 
different subbands lying in T and X valleys for a 33 atomic 
layer thick Si film. 



film (~ 1.086nm), only T± and T2 contribute, the spatial 
distribution of the total electron concentration is dictated 
by wavefunction distribution of only these two subands. 
The peak concentration comes at the middle of the film, 
and reduces as it approaches the surface. However, for 
a 33 monolayer thick film (~ 4.344nm), 4 T subbands 
and bottom most X subband contribute, and this in turn 
affects the total carrier distribution, as shown in Fig. [HJ 

3.3 Parabolic effective Mass and It's Va- 
lidity 

A simple parabolic effective mass has been derived in this 
section at the minima of different subbands to show some 
interesting transport properties of thin film. Parabolic 
effective mass m*(i,j) for the i th valley and j th subband 
is defined as 




[01] 



Figure 9: Variation of parabolic effective mass at the min- 
imum of first r subband with crystal direction and Si film 
thickness. Anisotropy is observed at small film thickness. 
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m*(i,j) 



ri 2 



d^E(i,j)/dk 2 (i,j) 



(4) 




Fig. [9] and 10 show how parabolic effective mass (m*), 
normalized to electron rest mass (mo), calculated at the 
minima of two bottom most subbands Ti and T2 vary 
along different crystal direction, for four different film 
thicknesses. One should note that, in both cases, for 9 
monolayer thick film, effective mass is highly anisotropic. 
For Ti, effective mass increases as one moves from [10] 



Figure 10: Variation of parabolic effective mass at the 
minimum of second T subband with crystal direction and 
Si film thickness. For very small thickness (~ lnm), < 
111 > effective mass is smaller than < 100 > effective 
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direction to [11] direction, whereas, it reduces for r 2 val- 
ley. However, for larger thickness, effective mass in both 
the valleys, becomes fairly isotropic. Since the effective 
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Figure 11: Distribution of < 100 > M*(z) along film 
thickness for various values of Si film thickness. Increas- 
ing thickness increases M* (z) and reduces spatial unifor- 
mity. 

masses vary with subbands, and electron concentration in 
different subbands has different spatial distributions, it is 
expected that effective mass should also have a spatial 
distribution. A 'distributed effective mass', say M*(z), a 
function of the depth z along the thickness of the film, 
has been defined as: 



M*{z) 



1 



E 



(5) 



where Wij(z) represents the fractional contribution of 
electron concentration at depth z from j th subband of 
i th valley. This way of defining < 100 > M* (z) has the 
underlying assumption that all the electrons (more gen- 
erally, an equal fraction of electrons from each subband 
of every valley) are moving along < 100 > direction. Fig. 
11 shows that for thinner films, < 100 > M*(z) is more or 
less uniform (which is because all the electrons are in T\ 
and r 2 subbands possessing almost same < 100 > effec- 
tive mass), but as thickness increases, the effective mass 
at position closer to surface becomes larger than that of 
the central part of the film. Thus, for larger thickness, 
carriers closer to center of the film are expected to be 
more mobile than those which are closer to the surface. 
Note that, this effect is inherent to the intrinsic film, com- 
ing from spatial distribution of wave functions associated 



with different subbands. Another interesting observation 
is that, a 9 monolayer thick film has (merginally) larger 
< 100 > effective mass than a 17 monolayer thick one at 
all z. This is because of the fact that the r 2 valley has 



larger effective mass for 9 monolayer thick film (Fig. 10 1 
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Figure 12: Original tight-binding bandstructure data and 
corresponding parabolic fits at four bottom most T sub- 
bands and bottom most X subband for < 100 > Si. The 
fits are reasonable for energy values less than ~ 0.5eV 
referenced from corresponding subband minima. 



However, one should note that, parabolic effective mass 
approximation is valid only for electrons with smaller en- 
ergy in the conduction band. The solid curves in Fig. [12] 
show the E — k relationship along the < 100 > direction 
for the four bottom most conduction subbands, calcu- 
lated from tight-binding method, as described in sec. [2] 
The dotted curves show the fitted parabolic bands with 
same effective masses, as calculated from eqn. @. AE S 
represents the electronic energy range in the s th subband 
between which the parabolic E — k tracks tight-binding 
E 
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it is clearly visible that 



k fairly well. From Fig. 
parabolic bands fail to track the actual bands for elec- 
tron energies in excess of ~ 0.5eV, referenced from corre- 
sponding band minimum, in all the cases. Nevertheless, 
the above simple analysis gives good qualitative insight 
about transport and mobility. 
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4 Self-consistent Solution of 
Poisson-Schrodinger Equation 

To study the effect of gate voltage on this structure, a 
model device has been assumed consisting of a top gate 
and a bottom gate separated from the film by thin insu- 
lator layers, as shown in Fig. [T] The Si film is assumed 
to be undoped. Thus channel charge corresponds to only 
mobile charge. If 4>{z) is the potential at z, then one can 
write the 1-D Poisson equation as [3J [TU] 



d 2 cj>{z) _ qN 
dz 2 eoe rs i 



e k B T (6) 



where q is electronic charge, e is permittivity of vacuum 
and e rs is the relative permittivity of the channel mate- 
rial. The boundary conditions are derived from the fact 
that the normal component of the displacement vectors 
inside Si and insulators will be the same at z = and 
z = t. Thus, at the boundaries, one can have 



and 



d<Kz) 
' dz 

d(j>(z) 



V B i - V, 



m 



e r2 



V g2 - V fb2 - 



(7) 



(8) 



dz " " t ox2 
e r \ and e r2 are the relative permittivities of insulatorl 
and insulator2 respectively, and, t ox \ and t ox2 are corre- 
sponding thickness of the insulators. Vfu is the flatband 
voltage between the i th gate and channel. Eqn. ^ can 
be solved iteratively to find <j)(z). To do this, first one 
assumes an initial potential profile 4>{z), and then calcu- 
lates bandstructure, which in turn provides the correction 
to (j>{z). This is iterated until it converges. However, one 
should note that, in every iteration, one needs to calculate 
bandstructure. This is because the potential <j>(z) adds a 
z dependent perturbation to the crystal potential, and 
thus both Ej and ip* are function of 4>(z). This makes 
the problem computationally intensive. 

To reduce computation, the following approximation 
has been made. Note that, the potential <f>(z) does not 
change drastically along z (which is shown later), and 
thus, as far as change in bandstructure is concerned, it's 
a fair assumption, that <j>(z) is constant (=<fi c ) along z. 
Suppose, H is the original unperturbed 107V x ION tight- 
binding Hamiltonian, E and ipo are the unperturbed 
eigen values and eigen functions respectively. If one as- 
sumes that external potential only changes the on-site 
energies, and not overlap integrals, then the perturbation 
AH can be written as 



Thickness 


terror (%) 


Tlerror (%) 


(AL) 


Mean 


SD 


Mean 


SD 


9 


-0.14 


0.03 


0.51 


2.09 


33 


0.65 


0.31 


2.98 


6.28 



Table 1: Mean and Standard deviation values of percent- 
age error in <j)(z) and n(z) for 9 and 33 monolayer thick 
films with V g = 1.5V. 



where / is 107V x lOiV diagonal unity matrix. Then, 

HiP = (H + AH)tJj = (E - q<f>c)ipo (10) 

which means that all the energy eigenvalues will be shifted 
by same energy (in other words, no relative change in 
energy eigenvalues) , and the wave functions remain in the 
unperturbed state. Thus, it is sufficient to calculate the 
bandstructure only once, and the same Ej and ipj can 
be used through all the iterations. This reduces the total 
runtime by nearly same number of times as the number 
of iterations it takes to solve the Poisson equation (which 
varies roughly from 20 to 200 for different cases). 

To validate the approximation, the amount of error be- 
ing incurred in the worst case (maximum gate voltage 
where band bending is maximum) has been examined and 
the results are tabulated in Table [TJ for both 9 and 33 
monolayer thick films. The error here has been defined as 
Perror = ^'""p""^"" x 100% where Papprox is the ap- 
proximate value of parameter P and P exa ct is the exact 
value of the parameter. 

To simplify the analysis, in the following, the metals 
used as gate electrodes, have been assumed to have mid- 
gap work-function, and charge trapping inside insulators 
is taken to be zero. This essentially means that the flat- 
band voltages Vfbi and Vfb 2 are taken to be zero. For 
simulation, it has been assumed that t ox \ — t ox2 = lnm, 
V g i = V g2 and e r i = e r2 = 3.9. 

4.1 Volume Inversion 

Under the above mentioned assumption that bandstruc- 
ture remains fairly the same under application of gate 
voltage, one can write the carrier density distribution n(z) 
as 



n(z) = n u (z)e k s T 



(11) 



AH = -q<b r I 



where n (z) is the intrinsic carrier density at z and is 
given by eqn. ([3| and <j>(z) is the potential profile ob- 
tained by solving the eqn. ([fj]). Fig. 13 and 14 show the 
(9) distribution of carrier density over different subbands and 
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the total density, for 9 and 33 atomic layer thick films re- 
spectively, when a IV supply has been applied to both the 
gates. One should note the difference in shape of the car- 
rier distribution as compared to intrinsic case. Also, the 9 
atomic layer thick film shows a higher peak carrier density 
compared to the 33 atomic layer thick one, although the 
total integrated carrier concentration is larger for thicker 
film. This can be explained with 'potential pinning' ef- 
fect, as discussed later. However, at any z, the fractional 
contribution from a subband to the total carrier density 
remains the same at any applied gate voltage. Fig. 15 




Figure 13: Total electron density distribution along the 
thickness and contribution from different subbands for a 
9 atomic layer thick Si film with V g — 1.0V. 



and 16 show the normalized carrier distribution along the 



depth z at different gate voltages, for 9 and 33 atomic 
layer thick films. The corresponding potential distribu- 
tion plot is shown in Fig. [17] As gate voltage increases, 
the carriers of the 'central-peaked ' channel start spreading 
toward the channel-insulator interface, and beyond a par- 
ticular threshold gate voltage V v t, the carrier distribution 
will start showing two 'humps' representing separation of 
peaks of carrier concentration inside channel. V vt is de- 
fined by the condition 



d 2 n{z 



dz 2 



v„=v vt 







(12) 



V v t for a thinner film is expected to be larger than a 
thicker one. In other words, carriers try to stay closer to 
the center for thinner film. With increase in gate volt- 
age, the shifting of carrier density peaks toward the gates 
can be thought of as a 'carrier pulling effect ' of the gate 




Figure 14: Total electron density distribution along the 
thickness and contribution from different subbands for a 
33 atomic layer thick Si film with V g = 1.0V. 



voltage. The effect can be explained from the potential 
profile <j>(z) which is smaller at points closer to the center 
of the film. From eqn. (11), one notices that the overall 



carrier density profile arises from the individual contri- 
ve*) 

bution of the two terms n (z) and e k B T . n (z) peaks 

at the center of the film, and reduces toward the surface 

»»(») 

whereas e k s T follows the opposite trend. Thus, at suffi- 
ciently large voltage, it is possible that the peak of carrier 
density occurs at the surface, which is qualitatively same 
as the classical picture, where the exponential term dom- 
inates so much that the effect of 'quantum mechanical' 
distribution of intrinsic carrier density gets nullified. As 
evident from Fig. |17| the potential profile for thinner 
films is more uniform, and actually 'pinned' at a higher 
value compared to the films of larger thickness leading 
to higher peak carrier concentration in thinner films as 
shown in Fig. [13] and [M] 

The above explanation becomes even more clarified 
from Fig. 18 Initially, for small gate voltage, both the 
surface potential <j) s and film center potential <po increase 
simultaneously at the same rate, and thus there will be 
only a single channel whose peak is at the center of the 
film. Beyond a certain gate voltage, <f> s and 4*o bifur- 
cate, and <po saturates very quickly. However, <fi s keeps 
increasing, though at a much slower rate than earlier, 
causing higher carrier concentration at points closer to 
the surface, finally destroying the single peaked channel 
and creating two channels of 'double hump' shape. Note 
that, for 9 atomic layer thick film, the 'pinning ' voltages 
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Figure 15: Normalized total electron density distribution 
along the film thickness for a 9 atomic layer thick Si film 
with V g — 0.5V, 1.0V and 1.5V. Single peak is observed 
even at larger gate voltages. 
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Figure 17: Potential distribution <f>(z) along the thickness 
of thin Si film with three different gate voltages: 0.5V, 
1.0V and 1.5V. The solid and dotted curves represent 33 
and 9 monolayer thick films, respectively. 



are higher than 33 atomic layer film. 




0.4 0.6 
Normalized z 



Figure 16: Normalized total electron density distribution 
along the film thickness for a 33 atomic layer thick Si film 
with V g = 0.5V, 1.0V and 1.5V. 'Double hump' charac- 
teristics is clear at higher gate voltages. 




Figure 18: Variation of surface potential <f> s (solid lines) 
and film center potential <f>o (dotted lines) with applied 
gate voltage V g , for two different film thicknesses. 



4.2 Channel Charge Distribution 

The 'carrier pulling effect ' explained in the previous sec- 
tion can have immense impact on the spatial distribution 
of total channel charge and hence device performance. 
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Consider a DGFET where the source end is grounded 
and the drain is connected to Vdd which is same as ap- 
plied gate voltage V g . At any point x along the channel 
(with x = being taken as the source end), suppose the 
quasi- Fermi level is V q f(x). V q f(x) can be assumed as 
independent of z. Then, the Poisson equation in eqn. Q 
gets modified as [5] 



d 2 4>(x, z) 
dz 2 



qN 
€Qe rs t 



EE 2 /(4)i#( 



J k 



xe 



(13) 



where all references have been made from grounded 
source chemical potential. The carrier density n(x, z) 
should now be a function of V q f{x) as well, which in turn 
depends on the drain voltage. Deriving the exact drain 
current for a nano-MOSFET needs proper attention on 
transport model and the details will be communicated 
separately. However, to get a quantitative estimate, a 
long channel device (L = 1/im, W = Ifim) with con- 
stant mobility has been assumed. The drain current and 
V q f(x) have been calculated using a similar procedure as 
in [6], by self-consistently solving eqn. (13 1 with drain 
current continuity equation. The extracted normalized 
carrier concentration n'{x,z) is plotted over the whole 
channel region in Fig. [19] for different cases. At any 
x = xq, n'(xo, z) has been defined as 



n'(x ,z) = 



n(x 0l z) 
MAX z {n{x ,z)} 



where MAX Z {.} represents the maximum value of {.} 
over z. The normalization is done in such a way which 
clearly shows the spatial shape of carrier distribution at 
different x. It is observed that to the source end (x = 0), 
due to larger potential difference between gate and chan- 
nel, there clearly exist two distinct carrier density peaks. 
However, as one moves toward the drain end, the po- 
tential difference between gate and channel reduces, and 
the two distinct peaks merge together producing a sin- 
gle center-peaked channel. Also, The magnitude of the 
total carrier concentration reduces toward the drain end. 
Putting in another way, near the source, carriers stay 
closer to the surface, and near the drain, carriers stay 
closer to the center of the channel. Thus, near the source, 
one expects more surface scattering and away from it, sur- 
face scattering is expected to reduce. Similar effects are 
true for gate leakage and gate capacitance, which can no 
longer be assumed to be uniform along the channel. It is 
evident from Fig. [19] that this effect is more pronounced 



in thicker channel devices, and at higher operating volt- 
ages. If the channel is thin enough, as is the case of (b2) 
in Fig. [19J it is possible to have a single peaked channel 
all over the device. 



(a1) 



(a2) 




Figure 19: Normalized carrier distribution over the whole 
Si channel for four different cases: (al) 33 monolayer thick 
channel, V g = 1.5V, Vdd = 1.5V (a2) 33 monolayer thick 
channel, V g = 1.0V, Vdd = 1.0V (bl) 9 monolayer thick 
channel, V g = 1.5V, Vdd — 1.5V (b2) 9 monolayer thick 
channel, V g = 1.0V, Vdd = 1-0V. 



( 14 ) 4.3 Evolution of Effective Mass Along 
DGFET Channel 



In sec. [3] it has been discussed in detail how < 100 > 
M* (z) varies along the thickness for a thin film Si. Now, 
when one considers the channel of a DGFET, as has been 
discussed in the previous section, the potential difference 
between gate and channel changes from source end to 
drain end. Thus total number of carriers at different sub- 
bands also changes along the channel. Keeping this in 
mind, one can define an 'average effective mass', M*(x) 



M* e {x) 



Jq n(x, z)dz 
Jo 



MHz) 



dz 



(15) 



which is basically harmonic average over the carriers 
along the thickness at a particular position x along chan- 
nel length. Physically, this indicates the 'average' effec- 
tive mass of an electron located at distance x from the 



source along the channel. Strictly speaking, eqn. ( 15 ) is 



valid only under the assumption that the vertical field is 
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fairly constant and each electron suffers same scattering where Q S i{z) is given by 
rate. Although this is not a very good approximation, but 
it gives an idea of how the channel charge distribution 
can affect spatial distribution of carrier effective mass. Qsi(z) 
Fig. 



20 shows for very small thickness channel (e.g. 9 



qN 



££2/(4)iv>f( 



(17) 



Fig. [21] shows the variations of channel capacitance with 



V =1.5V 



0.28 : 
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0.24 9 



0.22 
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Figure 20: Variation of 'average effective mass' M*[x) of 

carriers with channel along < 100 > direction for different Figure 21: C sl as a function of V g and depth z along 
channel thickness values and different gate voltages V g thickness of channel for (a) 9 monolayer thick channel 



with Vdd = Vg 



and (b) 33 monolayer thick channel. 



monolayer thick) , M* (x) hardly varies with x as well as 
V g . However, for higher thickness, a gradual decrease in 
M* (x) is observed from source end to drain end, and the 
effect is more prominent for higher gate voltages. 



z and V g . One should note that the position z = z max 
inside the film, where the coupling with gate is maximum, 
varies with gate voltage and as gate voltage increases, 
Zmax shifts toward the surface. 



4.4 Channel Capacitance 

Qualitatively, compared to classical analysis, the charge 
distribution in a quantum analysis, has two major dif- 
ferences: 1) The total charge in the channel reduces and 
2) The charge distribution peak shifts from the surface 
toward the center of the film. The extent of the shift 
depends on the applied gate voltage. Both these effects 
cause a change in total gate capacitance [5] . The channel 
capacitance per unit volume C' S i(z) at a depth z can be 
defined as the rate of change of charge per unit volume 
with respect to the potential at that point. Mathemati- 
cally, 



C si (z) 



dQsi(z) 
d<f>(z) 



(16) 



5 Conclusion 

A detailed analysis of generic ultra-thin-body DGFET 
has been performed in this work. The channel material 
has been chosen to be Si, but the analysis and methodol- 
ogy can be readily extended to other promising channel 
materials as well. Ultra thin film of Si has been shown 
to have larger and direct band gap, as opposed to bulk. 
It has also been shown that the intrinsic carrier concen- 
tration is not only less compared to bulk, but also has 
a distribution over the channel thickness, peaking at the 
center. The contributions of different subbands from dif- 
ferent valleys to both intrinsic carrier concentration as 
well as effective mass have been analyzed. The spatial 
distribution of distributed < 100 > effective mass along 
thickness has been studied. It has also been shown that 
along < 100 > direction, parabolic effective mass is fairly 
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valid till an electronic energy of ~ 0.5 eV from the corre- 
sponding subband minima in all the relevant valleys. A 
detailed insightful analysis of volume inversion has been 
performed and using 'carrier pulling effect ' of gate volt- 
age, channel charge distribution in a DGFET has been 
predicted. The effects of channel charge distribution on 
effective mass and channel capacitance have been ana- 
lyzed critically. 
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